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ABSTRACT 

Globular clusters (GCs) are the oldest stellar system in the Galaxy, in which the 
millisecond pulsars are widely believed to be the only steady 7 -ray emitters. So far 
9 7 -ray GCs have been identified and a few candidates such as 2 MS-GC 01 and IC 
1257 have been suggested. In this work, after analyzing the publicly-available Fermi 
LAT data we confirm the significant 7 -ray emission from 2MS-GC01 and IC 1257 and 
report the discovery for 7 -ray emission from NGC 5904 and NGC 6656 within their 
tidal radii. Also a strong evidence of significant 7 -ray emission is found from FSR 
1735. From the observed 7 -ray luminosities, the numbers of MSPs that are expected 
to be present in these GCs are estimated. 

Key words: Gamma rays: observations — Globular clusters: individual (NGC 5904, 
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1 INTRODUCTION 

Globular clusters (GCs) are spherical stellar systems in 
which stars are tightly bound by gravity. They are the oldest 
and also the most dense stellar systems in the Galaxy. On 
average the stellar density is ~ 0.4 star per cubic parsec, in¬ 
creasing to 100 or 1000 stars per cubic parsec in their cores. 
The dynamical interactions with high stellar encounter rate 
are expected to form low-mass X-ray binaries (LMXBs). The 
formation rate per unit mass of LMXBs in GCs is known to 
be greater t han in other places of the Galaxy (|Kate |[l975l : 
IClark |[l975h and millisecond pulsars (MSPs) are generally 
belie ved to be the descendants of LMXBs ([Alpar et al. I 
Il982h . It is also widely accepted that MSPs contribute al¬ 
most all 7 -ray emissions of GCs. Before 2008 no GCs had 
been known as 7 -ray emitter. The situation changed dramat¬ 
ically after the successf ul launch of Fermi Gamma-ray Space 
Telescope (| Atwood et al. 1 120091 since the Large Area Tele¬ 
scope (LAT) onboard the Fermi satellite, a pair-production 
telescope, can detect 7 -rays at energies between ~ 20 MeV 
and > 300 GeV with unprecedent sensitivity. The discovery 
of 7 -ray emission from GCs was made in 47 Tucanae (NGC 
104; Abdo et al. 2009), in which at least 23 MSPs have 
been known by radio/X-ray observations. As suggested in 
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lAbdo et al. I ( 20101 ), the 7 -ray emission from GCs is not due 
to single MSP but all the MSPs as a population. 


Mor e than 150 GCs have been identified in optical band 
([Harris I Il996h , but most of them are not significant 7 - 
ray emitters. Nine 2FGL sour ces are do ubt less a ssociated 
with GCs, including 47 Tuc (lAbdo et al. 2009 k Omega 
Cen (NGC 5139), NGC 6266 NGC 6388 Terzan 5, NGC 
6440, NGC 6626, NGC 6652 (lAbdo et al. Il 2010 h and M 80 
(NGC 6093; iTam et al. 1 12011 ). Moreover, th ere are a few 
cand idates such as 2MS-GC01 and IC 1257 ([Nolan et al. I 
l2012h . The non-detection of significant 7 -ray emission from 
most GCs is likely due to their low luminosities of the GeV 
emission and due to the strong background emission for the 
global clusters superposed in the Galactic disk. 


In this work, we analyze the latest Fermi LAT data 
aiming to identify more 7 -ray emitters of GCs. This work is 
structured as the following. In section 2 we present our data 
analysis. Our results are introduced in Section 3. Section 4 
is our summary with some discussion. 
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2 DATA ANALYSIS 

2.1 Data preparation 

We use more than 6 years of Fermi LAT data (Pass 7 Re¬ 
processed data; totally covering 2191 days from August 4th 
2008) with energy ranges from 100 MeV to 100 GeV. The 
Fermi Science Tools v9r33p0 package is used to analyze the 
data, which is available from the Fermi Science Support Cen- 
teiB The P7REP.SOURCE (evclass=2) photons for point source 
analysis are selected in this work. To reduce the contribution 
from Earth albedo 7 -rays, we only use events with zenith 
angles ^ 100°. The P7REP_S0URCE_V15 instrument response 
functions (IRFs) is used. 

2.2 Analysis method 

We try an overall search of the GCs listed in the catalogue 
compiled by Harries (1996). In the analysis of each GC, pho¬ 
tons from a 20° x 20° square region of interest (ROI) are 
selected, and binned into spatial pixels of 0.1° X 0.1°. The 
first step is to apply standard binned maximum likelihood 
analysis using gtlike. We modelled the events with compo¬ 
nents of target source (with coordinates given in the cat¬ 
alogue) and backgrounds. The background components in 
this model are composed of all sources in the second Fermi 
LAT catalogu43 within the ROI and diffuse components, in¬ 
cluding the Galactic diffuse model (gll_iem_v05_revl. f it) 
and isotropic background (iso_source_v05_revl.txt). As 
usual, a single power law (PL) is used to model the point 
source. We also apply an exponential cutoff power law (PLE) 
for the 7 -ray GC candidates, in which radio pulsars have 
been detected. For both models, the prefactor and the spec¬ 
tral index are set free. Additionally, the cut off energy is set 
free in PLE model. 

We then choose the GCs with Test Statistic (TS) values 
greater than 25 for a further study. The TS is defined as TS 
= 2(£1 — £ 0 ), where £0 is the logarithmic maximum likeli¬ 
hood value for a model without the source (null hypothesis) 
and £1 for a model with additional source at a specified 
location. For those 7 -ray GC candidates with TS ^ 25, we 
created 5° x 5° residual TS maps for them using the tool 
gttsmap. The TS maps aiming to identify weaker sources are 
created by moving a putative point source through a grid of 
locations on the sky and maximizing -log(Likelihood) at each 
grid poin 1 @. In taking this step, the target source that cor¬ 
responds to the GC is removed from the model. All point 
source parameters are fixed while the diffuse components 
are set free. To reduce the contamination from background 
emissions due to point spread function (PSF) at lower en¬ 
ergy band, different energy ranges for each GC are selected 
in TS map creation. 

From the residual TS maps, we find NGC 5904, NGC 
6656 have evidences of 7 -ray emission within their tidal 
radii. Meanwhile, there is a 7 -ray emission evidence from 
FSR 1735 (see lower panel in Fig. G2>, the separation be¬ 
tween the best-fit position of emission and core of FSR 

1 http://fer mi. gsfc.nasa.gov/ssc/data/analysis/software/ 

2 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/ 

3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/likelihood_t 

utorial.html 


1735 is 15.6'. Furthermore, iNolan et al. I (l2012h reported 
other 2 associations with 2FGL sources in second cata¬ 
logue: IC1257 (2FGL J1727.1-0704) and 2MS-GC01 (2FGL 
J1808.6-1950c). We also create TS maps for them. We then 
use gtfindsrc to search for the positions of these 7 -ray can¬ 
didates. 

The significance from corresponding TS values follows 
the x 1 2 distribution with 2 degrees of freedom (PL model; 3 
degrees for PLE model). To avoid fake signal with a num¬ 
ber of se arch positions, sim ilar estimation of trial factor is 
involved (|Tam et al. I [ 2011 ). Totally 35 GCs with TS > 25 
are further analysed in this work, so Ngc — 35. The num¬ 
ber of bins, Nbin ~ 15, is derived by dividing the averaged 
tidal radius area of all known GCs into 0.1° x 0.1° pixel, 
where the average tidal radius is 13k We get the post-trial 
significance using (1 — pY = 1 — p', where p is the pre-trial 
p- value, / = Ngc x Nbin is the trial factor, and p is the 
post-trial p- value. 

Finally, we obtain the spectral energy distributions 
(SED) of these GCs. We divide the energy range from 200 
MeV (100 MeV for 2MS-GC01) to 100 GeV into logarithmi¬ 
cally equally spaced energy bins (base = 1.94). The flux is 
obtained by fitting all model components in each bin using 
gtlike. In the source model, the normalizations of Galactic 
diffuse component and isotropic component are left free. 

3 RESULTS 

3.1 2MS-GC01 

2MS-GC01 h as formally been associated with 2 FGL 
J1808.6-1950C llNolan et, al. 1120121: ICholis et, al. 112014l . Our 
best-fit position is R.A. = 272.18, decl. = -19.93 (J2000), 
which is offset by 4.8 ; from the GC core (see the TS map, 
left panel of Fig. We obtain an index of 2.2 db 0.1 
and a cutoff energy E c =3.4 =L 0.2 GeV in a PLE spec¬ 
tral fitting. A TS value of 241, corresponding to a detec¬ 
tion significance of 15.1cr (14.8a post-trial) is given. If we 
just take into account the photons above 400 MeV (energy 
range selection for creating TS map), the detection signif¬ 
icance is ~ 10 a. Fig. [T| shows the spectrum of the 7 -ray 
emission in E 2 * dN/dE. The integrated photon flux in the 
energy of 100 MeV—100 GeV is Fo.i-iooGeV = (12.4 ± 
1 . 2 ) x 10 8 cm 2 s 1 and the integrated energy flux is 
Fo.i-iooGeV = (5 .7 =fc 0 .5 ) x 1 0 ~ n erg cm -2 s -1 . At a dis¬ 
tance of 3.6 kpc ([Harris ll 19961 . 2010), the 7 -ray luminosity 
is Lo.i-iooGeV = (9.2 db 0.8) x 10 34 erg s -1 . 

3.2 IC 1257 

IC 12 57 has formally be e n associated with 2 FGL J1727.1- 
0704 (jNolan et al. II 2 OI 2 I ; ICholis et al. 1 120141 5 . Our best-fit 
position is R.A. = 261.98, decl. = -7.07 (J2000), which is 
offset by 10 . 6 / from the cluster core (see the right panel in 
Fig. [2]). Fitting the spectrum with a PLE model we have 
an index of 1.2 ± 0.7 with a cutoff energy E c = 1.6 ± 0.9 
GeV. We obtained a TS value of 35 and the corresponding 
detection significance of 5 . 2 cr (4.0cr post-trial). We have TS 
= 30 if instead just the photons above 400 MeV have been 
taken into account, corresponding to a significance of 4.7cr 
(3.4cr post-trial). It is also the energy range selection for TS 
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map creating. The integrated photon flux between 100 MeV 
and 100 GeV is Fo.i-iooGev = (5.7 zb 1.0) x 10 -8 cm _ 2 s _1 , 
while the integrated energy flux is Fo.i-iooGeV = (5.0 zb 
0.9) x 10 -12 erg cm -2 s -1 . The luminosity is Lo.i-iooGeV — 
(3.7 =b 0.7) x 10 35 erg s _1 at a distance of 25 kpc. Such a 
luminosity is rather high in comparison with other 7 -ray 
GCs, which is actually the largest among all detected 7 -ray 
emitting GCs by now. 

3.3 NGC 5904 

In this GC five radio pulsarfl have been detected 
([Anderson et al. I Il997l : iHessels et al. I I2007T ) but no signif¬ 
icant 7 -ray emission has been reported. In our analysis we 
find some evidence for 7 -ray emission from NGC 5904. With 
the photons in the energy range of 0.1 — 100 GeV, the TS 
value under PL model is 26, corresponding to a significance 
of 4.6cr (see the upper left panel in Fig. [3|) and the post¬ 
trial significance is 3.2cr. To minimize the influence of strong 
background emission at low energies, we choose the photons 
at energies of ^ 500 MeV (for these photons the PSF is much 
smaller than that at energies of ~ 100 MeV and the back¬ 
ground pollution can be effectively suppressed) to derive TS 
map, which gives a TS value of 21 (~4.0a; post-trial 2.4a). 
Our best-fit position is R.A. = 229.69, decl. = 2.17 (J2000), 
5.9' from the core position of N GC 59 04, well within the 
tidal radius of NGC 5904 ([Harris 1996 . the 2003 version). 
The SED can be well fit by a single power law with an index 
of 2.3 =b 0.2. The integrated photon flux is Fo.i-iooGeV = 
(6.2 zb 1.4) x 10 -9 cm -2 s -1 and the integrated energy flux 
is Eo.i-iooGeV = (3.5 ± 0.8) x 10 -12 erg cm -2 s _1 . When 
deriving TS map, an additional source (R.A. = 230.54, decl. 
css 4.32) is added to the model, of which the spectrum is set 
as a single power law and to the ROI center with a sepa¬ 
ration of 2.4°. At a distance of 7.5 kpc, the 7 -ray emission 
luminosity is Lo.i-iooGeV = (2.4 zb 0.5) x 10 34 erg s -1 . We 
also use a PLE model to fit the spectrum, where the photon 
index F = 1.2 zb 1.3 and the cutoff energy is (2.6 ± 2.8) GeV. 
In this model, the significance keeps on the same level, 3.1a 
post-trial. No significance of the PLE model over the PL 
model, so we consider the PL model as a better model. 

3.4 FSR 1735 

The 7 -ray emission from FSR 1735 is found with a TS value 
of 55 and the corresponding significance is ~ 7.0a (6.2a 
post-trial). Again, we only take into account the photons 
above 400 MeV to create TS map (see lower panel in Fig. [3j) , 
the value of TS is 39 and the corresponding significance is 
5.8a (4.8a post-trial). Our best-fit position is R.A. = 253.37, 
decl. = —46.91 (J2000). Though the tidal radius of FSR 
1735 is unknown, the bset-fit 7 -ray position found in our 
analysis just has an offset ~ 15.6 X from the core of FSR 1735. 
Using a single power law model, we obtained a photon index 
T — 2.2 zb 0.1. The integrated photon flux is Fo.i-iooGeV = 
(2.7 zb 0.4) x 10 -8 cm _ 2 s _1 and the integrated energy flux 
is Fo.i-iooGeV = (1.8 zb 0.3) x 10 -11 erg cm -2 s -1 . At a 
distance of 9.8 kpc, the 7 -ray luminosity is Lo.i-iooGeV = 
(2.1 zb 0.3) x 10 35 erg s" 1 . 

4 http://www.naic.edu/ pfreire/GCpsr.html 


3.5 NGC 6656 

Two radio pulsar^ have been discovered ([Lynch et al. I 
B) in this GC. It also has no association in the second 
catalogue of LAT. Our reduction on this region show an evi¬ 
dence of 7 -ray emission. The best-fit position (R.A.=279.30 
decl.=-24.05; J2000) is well within the tidal radius circle of 
NGC 6656. Using a single power law model (3.4a post-trial), 
we obtained a photon index of T = 2.7 zb 0.1, and 0.1-100 
GeV photon and energy fluxes of (2.3 zb 0.6) X 10 -8 cm -2 
s -1 and ( 8.6 zb 1.9) x 10 “ 12 erg cm -2 s -1 , respectively. The 
7 -ray luminosity is Lo.i-iooGeV — (1.1 zb 0.2) x 10 34 erg s -1 
at a distance of 3.2 kpc. The spectrum also can be fitted with 
a PLE model, where the photon index T = 2.0 zb 0.5 and 
cutoff energy E c = (2.4 zb 1.8) GeV. Due to the significance 
under PLE model is also 3.4a post-trial, so we consider PL 
model as the better one. 


4 SUMMARY 

In this work we have analyzed the Fermi LAT data of more 
than 150 GCs and detect or find some evidence of 7 -ray 
emission from 3 GCs, i.e. FSR 1735, NGC 5904 and NGC 
6656. Another two (2MS-GC01 and IC 1257) have been re - 
ported as 7 -ray emission candidates in I Nolan et al. (|2012h . 
of which the associations of significant 7 -ray emission have 
been confirmed here. Among the three new detected can¬ 
didates, the evidence of 7 -ray emission from FSR 1735 is 
very significant (the post-trial significance is at a confidence 
level of > 5a). Though the tidal radius of such a GC is un¬ 
known, the 7 -ray source found in our analysis just has an 
offset ~ 15. 6' from the core of FSR 1735, which is small 
enough to allow us to suggest an association with these two 
sources. The 7 -ray emission signals from the other two GCs 
are weaker (the confidence level is above 3a but smaller than 
5a) and more data is needed to firmly establish their 7 -ray 
emitter nature. Fermi satellite continues to collect the 7 -ray 
data and more 7 -ray GCs are expected to be identified. Nev¬ 
ertheless for the GCs within the Galactic plane the prospect 
is not very promising due to the strong background that can 
minimize the 7 -ray signal effectively. More GCs may be de¬ 
tected in the future with a high angular resolution telescope 
dedicated to the sub-TeV (from ~ 10 MeV to ~ 10 GeV) 
7 -ray photon detection. For example, the proposed mission 
“PAir-productioN Gamma-ray Unit” aims to have an an¬ 
gular resolution a factor of ^ 5 better than the currently 
operating Fermi Gamma-ray Space Telescope in the sub- 
GeV range (Wu et al. 2014). The high quality 7 -ray maps 
will significantly improve the identification of “point-like” 
sources from extended and complicated diffuse 7 -ray back¬ 
ground, which is thus very suitable to identify more 7 -ray 
GCs. 

As already mentioned in Section 3, no pulsars have been 
reliably identified in these GCs except NGC 5904 and NGC 
6656. The detection of 7 -ray emission in turn provides us the 
chance to estimate the number of millisecond pulsars that 
are believed to be the unique stable 7 -ray emitters in GCs. 
Our approach is following Abdo et al. (2010) and the number 
of MSPs in each GC is estimated using their Eq.(l). Where 

5 http://www.naic.edu/ pfreire/GCpsr.html 
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(E) = (1.8 =b 0.7) xlO 34 erg s -1 is adopted as average spin- 
down power of pulsar, and the 7 -ray luminosity conversion 
efficiency is set to 0.08. The results are shown in Table [ 2 ] 
For NGC 5904 and NGC 6656, the estimated numbers of 
MSPs are more than that found currently, implying that in 
the future more powerful radio telescopes (or arrays) have 
very promising detection prospects as long as these MSPs 
also radiate non-ignorable amount of energy in radio bands. 
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Table 1. 


GC 

Name 

R.A. 

Position^ 1 ) 

Decl. 1. 

b. 

TidalC 2 ) 

radms(') 

Position^ 3 ) 

R.A. Decl. 

Offset 

(') 

NGC 5904 

229.64 

2.08 

3.86 

46.80 

28.40 

229.69 

2.17 

5.9 

NGC 6656 

279.10 

-23.90 

9.89 

-7.55 

28.97 

279.30 

-24.05 

14.2 

FSR 1735 

253.04 

-47.06 

339.18 

-1.85 


253.37 

-46.91 

15.6 

2 MS-GC 01 

272.09 

-19.83 

10.47 

0.10 


272.18 

-19.93 

4.8 

IC 1257 

261.79 

-7.09 

16.53 

15.14 


261.98 

-7.07 

10.6 


Notes. (1) Coordinates derived from catalogue of lHarrisl ( 19961 , 2003 version) and corresponding Galactic coordinates, ( 2 ) Tidal 
radius of FSR 1735, 2 MS-GC 01 and IC 1257 are not given, (3) Our best-fit position of 7 -ray emission from GCs. 


Table 2. Parameters of 7 -ray from GCs. 


GC 

Name 

Spectral 

model 

TS 

Significance 

(post-trial;cr) 

Photon 

index 

Cutoff 

(GeV) 

photon^ 1 ) 

flux 

Energy^ 2 ) 

flux 

Distance^ 3 ) 

(kpc) 

L( 4 ) 

N MSP 

NGC 5904 

PL 

26 

3.2 

2.3 ± 0.2 


6.2 ± 1.4 

3.5 ± 0.8 

7.5 

2.4 ± 0.5 

16 ± 4 


PLE 

28 

3.1 

1.2 ± 1.3 

2.6 ± 2.8 

1.6 ± 0.4 

1.9 ± 0.4 

7.5 

1.3 ± 0.3 

9 ± 2 

NGC 6656 

PL 

27 

3.4 

2.7 ± 0.1 


23.2 ± 5.9 

8.6 ± 1.9 

3.2 

1.1 ± 0.2 

7 ± 2 


PLE 

30 

3.4 

2.0 ± 0.5 

2.4 ± 1.8 

13.1 ± 4.1 

6.3 ± 1.4 

3.2 

0.8 ± 0.2 

5 ± 1 

FSR 1735 

PL 

55 

6.2 

2.2 ± 0.1 


26.9 d= 4.3 

18.4 ± 2.9 

9.8 

21.1 ± 3.3 

147 ± 23 

2 MS-GC01 

PLE 

241 

14.8 

2.2 ± 0.1 

3.4 d= 0.2 

124.7 ± 11.8 

56.8 ± 5.0 

3.6 

8.8 ± 0.8 

61 ± 5 

IC 1257 

PLE 

35 

4.0 

1.2 ± 0.7 

1.6 ± 0.9 

5.7 ± 1.0 

5.0 ± 0.9 

25.0 

37.4 ± 6.7 

260 ± 47 


Notes. ( 1 ) Integrated 0.1-100 GeV photon flux in unit of 10 9 cm 2 s 1 , ( 2 ) Integrated 0.1-100 GeV energy flux in unit of 10 12 erg 
cm -2 s -1 , (3) Distance from [Harris ( il996 . 2003 version), (4) 0 . 1-100 GeV luminosity in unit of 10 34 erg s 1 . 
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Figure 2. 5° X 5° TS maps of the ROIs for 2MS-GC01 and IC 1257. Both tidal radii are unknown. The cross represent the best-fit 
centriod of the 7 -ray emission from them. 



232.000 230.000 228.000 


TS value Right ascension 



0 5 10 15 20 



280.000 278.000 

Right ascension 



FSR 1735 J 

-45.000 


-46.000 


c 

o 

* 

1 -47.000 


Q 


-48.000 


-49.000 



256.000 254.000 252.000 250.000 


TS value 


0 5 


Right ascension 


10 15 20 25 30 


35 


Figure 3. 5° x 5° TS maps of the ROIs for NGC 5904, NGC 6656 and FSR 1735. The cross represent the best-fit centroid of the 7 -ray 
emission from them. The circles represent the tidal radii ( [Harris 19961 . 2003 version). The tidal radius of FSR 1735 is unknown. 
























